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The PHENIX experiment at RHIC measured u), r\ and X^-meson production at high px 
in p + p, d + Au and Au + Au collisions at y/sjviv = 200 GeV. Measurements performed 
in different hadronic decay channels give consistent results. The measured ratios of all 
three mesons to it are found to be fiat as a function of pt in p + p collisions and equal 
to ui/tt° = 0.81 ±0.02 ±0.07, 77/71-° = 0.48 ± 0.02 ± 0.02 and K° s /tt° = 0.45 ± 0.01 ± 0.05. 
Nuclear modification factor measured for oj-mesons in central Au+Au collisions is Raa = 
0.4 ± 0.15. 



1. Introduction 

Hadron spectra measured at high transverse momentum (px > 2 GeV/ c) pro- 
vide important information on particle production mechanism and properties of 
the matter produced in RHIC collisions. High-py hadrons are mainly produced by 
fragmentation of partons originating from hard scattering processes. In nucleon- 
nucleon collisions such processes are relatively well understood and corresponding 
cross sections can be calculated in perturbative QCD. In nucleus-nucleus collisions 
various nuclear effects like multiple scattering in the initial state, modification of 
the parton distribution and fragmentation functions, energy loss in the medium can 
modify particle production. These effects can be studied by comparison of particle 
properties measured in p + p and nucleus-nucleus collisions at the same energy. 
Measurement of u>, r\ and ii"g-meson production in p + p, d + Au and Au + Au 
collisions is a part of systematic study of particle properties at RHIC. 



2. Experimental setup and data samples 

Two central spectrometers of the PHENIX ^ experiment (Acf> < 90°, | i] |< 0.35) 
are used to measure neutral and charged particles produced in heavy ion collisions 
at RHIC. Beam-Beam Counters and Zero Degree Calorimeters provide the min- 
imum bias (MB) trigger, determine z-coordinate of the collision vertex and the 
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event centrality. Drift Chamber and the first layer of the Pad Chamber are used 
for reconstruction of charged particle momentum with resolution of <j(pt)/pt ~ 
1.0%pt © 1.1%. Two other layers of the Pad Chamber can be used for track confir- 
mation by matching it to the reconstructed hit. High segmentation electromagnetic 
calorimeter (EMCal) with energy resolution of a(E)/E sa 8.1%/s/E® 2.1% is used 
as a primary detector for reconstruction of 7-quanta and 7r°-mesons 

The experimental data samples gathered by PHENIX in years 2003-2005 and 
decays analyzed in these samples are summarized in Table[TJ Besides the MB trigger 
we use high-px 7-trigger (7) realized by adding together amplitudes in 4x 4 adjacent 
EMCal towers and comparing them to a threshold of 1.4 GeV in p + p and 2.4 GeV 
in d + Au collisions. Only events with the vertex |z„ e rt|< 30 cm from the center of 
the PHENIX experimental region and satisfying various quality assurance cuts are 
accepted in the analysis. 

Table 1. Analyzed decays and data samples. 

Run System Trig. Events /L Decays BR (%) p T (GeV/c) 

Run3 p + p 7 4.6 X 10 7 0.22 pb' 1 r) -> n°n+n- 22.7 ±0.4 3-8 

Lu^n°n+n- 89.1 ± 0.7 2.75- 9.25 

u> -* 7T°7 8.9±o.23 2.5 - 6.5 

K° s -> 7r°7r° 30.7 ±0.1 2.5-6.5 
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3. Analysis 

As the first step we reconstruct tt° — > 77 decay by combining EMCal clusters 
assuming that they are produced by photons originating from the collision vertex. 
To reduce number of fake 7r° candidates we select clusters with energy E n > 0.2 GeV 
in p + p and d + Au and £L > 0.35 GeV in Au + Au collisions. Shower profile cut is 
used to reject broader showers having mostly hadronic origin. In Au + Au collisions 
we use an asymmetry cut a =| _Ei 7 — E 2l | /(£a 7 + E 2l ) < 0.8 to further reduce 
background since high-py combinatorial pairs are peaked at a = 1 due to the 
steeply falling spectrum of single photons. Position and width of the reconstructed 
7r°-peak are parameterized as a function of px- The width of the 7r°-peak decreases 
from 12 MeV at p T = 1 GeV to 8 MeV at p T > 3 GeV/c. Pairs of photons 
with invariant mass within two standard deviations of the 7r°-peak position and 
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transverse momentum pr > 1 GeV/c are selected as 7r°-candidates. The candidates 
are assigned the exact mass of the 7r°-meson and measured pr of the pair. 

Depending on the decay mode we further combine a 7r°-candidates with other 
particles in the event. In case of u>(t)) — > 7r°7r + 7r _ decay we combine ir° candidates 
with any pair of negatively and positively charged tracks assuming them to be 
7r-mesons. Momenta of the charged tracks are selected to be below 4 GeV/c and 
8 GeV/c for 77 and w respectively. At higher momentum the charged particle spectra 
are dominated by decay products of long living particles. More details are described 
in SI. For uj — ► 7r°7 decay reconstructed 7r°-candidates are combined with all other 



photons from the same event with energy E 1 > 1 GeV. In case of K s 



decay 



we combine two 7r°-candidates together. However, if the invariant mass of two most 
energetic photons assigned to two different 7r°-candidates is within four standard 
deviations of 7r°-mass then such pair is rejected. This cut suppresses combinatorial 
background from misidentified 7r°'s. To account for K s -meson life time (cr = 2.7 
cm) the decay point for every pair is shifted from the collision vertex by the de- 
cay length of K s along the momentum vector of the pair. Finally we recalculate 
momenta of the 7r°-candidates and eventually momentum of the pair. 

Examples of invariant mass distributions measured in p + p collisions are shown 
in Figure [T] Mixed event technique does not reproduce background shape for none 
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Fig. 1. Invariant mass spectra for different decays at 5 < pxiGeV/c) < 6 in p + p collisions. 



of the decays because of particle correlations present in the events. To extract the 
particle raw yields we use several fitting options. First we fit the peaks with Gaussian 
function and adjacent background with a parabola as shown in Figure [U Secondly, 
the background is approximated in the same region with a parabola, but the points 
under the peak are rejected from the fit. The parabola interpolated under the peak 
is subtracted from the histogram. The yield is counted as a sum of the bins in the 
region initially excluded from the fit. Finally, we narrow the fitting region and use 
linear function instead of parabolic to describe the background. The three fits are 
applied to invariant mass spectra collected with different analysis cuts (matching 
of the charged track to the PC3 hit, tighter cuts for neutral particles) where the 
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shape of the background distribution is different. The raw yields extracted in each 
case are fully corrected for reconstruction and trigger efficiencies described further. 
The hrst of the described fits with its statistical error is taken as the result and a 
variance of the six measurements is added to the systematic errors as an estimate 
of yield extraction uncertainty. 

For the evaluation of meson reconstruction efficiencies we use Exodus as a par- 
ticle generator and PISA (PHENIX Integrated Simulation Application) for full 
Monte-Carlo simulation of the detector and online trigger settings. For the three 
body decay of u> and 77- mesons the EXODUS is also u sed as a particle decay engine 
weighted to reproduce the phase space density EEE3 The simulation code is tuned 
to represent the real configuration of the detector subsystems and to make sure 
that simulated shape, position and width of the 7r°, 77, uj and iTji-peaks are consis- 
tent with the ones measured in real data at all pr- Left panel of Figure [2] shows 




p T (GeV/c) p (GeV/c) 



Fig. 2. Geometrical acceptances and meson trigger efficiencies for hadron decays of u), Kg and 
r;-mcsons. Systematic uncertainties are shown with bands. 

the acceptance efficiencies which depend on the detector geometry, decay kinemat- 
ics, performance of detector subsystems and on the cuts used in the analysis. Right 
panel of Figure [2] shows the high-p T 7-trigger efficiencies for reconstructed decays in 
p+p. They are recalculated from the measured single photon trigger efficiency using 
kinematics of the particular decays. Both efficiencies must be taken into account 
to correct the raw yields measured in high-p^ 7-triggered data samples whereas 
only the curves shown in the left panel are used in MB samples. In Au + Au events 
one must account for the additional loss in efficiency caused by cluster overlap 
in the EMCal. This correction is based on simulated single w-mesons embedded 
in real events and reconstructed using the analysis procedure. Percentage of the 
reconstructed mesons is used to correct for this effect. 

Systematic errors for different decay modes of u> and K ^-mesons are summarized 
in Table [2 Details for r/-meson are described in 121 Total systematic error of the 
measurements is dominated by raw yield extraction, meson trigger efficiency and 
total cross section uncertainties. 
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Table 2. Relative systematic errors (%) for different decay modes and collision systems. Values 
with a range indicate variation of the systematic error over the px range of the measurement. 

Source lq — » tt°tt + it~ lu — » 7T°7 Kg — > 7r 7r" 

p + p d + Au p + P d + Au Au + Au P + P d + Au 

Acceptance 5 - 10 9- 12 10-20 8- 12 14- 16 10- 25 10 - 20 

Trigger efficiency 3- 10 5-7 2-7 5 - 2- 10 5 

Yield extraction 5 - 25 10- 15 5 - 15 10 15 - 35 7- 30 9 

MB trigger 10 8 10 8 4 10 8 

Total 15 - 25 18 - 22 15 - 25 17 - 20 20 - 45 20 - 40 15 - 25 



4. Results 



The invariant pt spectra measured for w-mesons at ^/snn = 200 GeV are shown 
in the left panel of Figure [3J In p + p collisions we observe a very good agreement 
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Fig. 3. Left: Invariant pt spectra for w-mesons in p + p, d + Au and Au + Au collisions at 
y/sjfN = 200 GeV. Right: (oj/tt°) ratio as a function of transverse momentum. Statistical errors 
are shown with error bars and systematic errors are shown with boxes. 



between Run3 and Run5 results and between the two decay modes measured 
in the analysis. In Run5 both the high-pr and the MB trigger data give consis- 
tent results. Dashed curves in d + Au and Au + Au points are the fits to p + p 
spectrum scaled by the corresponding number of binary collisions. As one can see 
the production of cj-mesons in d + Au and peripheral Au + Au collisions follows 
the binary scaling, while in the MB and in the most central Au + Au collisions 
we observe a suppression of w-meson production with Raa = 0.4 ± 0.15. Ratio of 
w-meson and 7r°-meson spectra ^ in p + p collisions is shown in the right panel of 
Figure [3l The ratio is flat in the pr range of the measurement. Fit to a constant 
shown with a solid line gives a value of 0.81 ± 0.02 ± 0.07 that agrees with previous 
PHENIX measurements of 0.85 ± 0.05 ± 0.09^. Both measurements give a value of 
the ratio slightly below PYTHIA ™ prediction of one. Lower energy measurements 
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of R-806121 and E706 ^ collaborations are consistent with the presented results. 

In Run3 we also measured 77 — > 7r°7r + 7r~ decay mode. The result is published 
in El. From the analysis standpoint this decay is very similar to the w-meson decay 
and the result is consistent with the 77 — > 77 measured by PHENIX. 

Invariant cross sections of the i^-meson production in p + p and d + Au colli- 
sions at y/SNN — 200 GeV are shown in the left panel of Figure 2] The K s — > 7r°7r° 
results are c omp leted with results of charged kaon measurements using PHENIX 
time-of-flightlSl. The Pt range of Run5 measurements extends to 10 GeV/ c provid- 
ing a baseline for heavy ion measurements. Cross section measured in p+p collisions 
is consistent with K s — ► 7r + 7r~ results from STAR in the overlap region 031. De- 
pendence of the K/tt ratio on transverse momentum is shown in the right panel of 
Figure |U Fit to a constant at p T > 2 GeV/c gives a value of 0.45 ± 0.01 ± 0.05. 
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Fig. 4. Left: Invariant cross section for it-mesons in p+p and d+Au collisions at y/Sf^N = 200 GeV . 
Right: K/ir° ratio as a function of transverse momentum. Statistical errors are shown with error 
bars and systematic errors are shown with boxes. 
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2. Major Headings 

Major headings should be typeset in boldface with the first letter of important 
words capitalized. 

2.1. Sub-headings 

Sub-headings should be typeset in boldface italic and capitalize the first letter of 
the first word only. Section number to be in boldface roman. 

2.1.1. Sub-subheadings 

Typeset sub-subheadings in medium face italic and capitalize the first letter of the 
first word only. Section numbers to be in roman. 

2.2. Numbering and spacing 

Sections, sub-sections and sub-subsections are numbered in Arabic. Use double 
spacing before all section headings, and single spacing after section headings. Flush 
left all paragraphs that follow after section headings. 

2.3. Lists of items 

Lists may be laid out with each item marked by a dot: 

• item one, 

• item two. 

Items may also be numbered in lowercase roman numerals: 

(i) item one 

(ii) item two 

(a) Lists within lists can be numbered with lowercase roman letters, 

(b) second item. 

3. Equations 

Displayed equations should be numbered consecutively in each section, with the 
number set flush right and enclosed in parentheses 



Equations should be referred to in abbreviated form, e.g. "Eq. |T|)" or "(2)". In 
multiple-line equations, the number should be given on the last line. 
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Displayed equations are to be centered on the page width. Standard English 
letters like x are to appear as x (italicized) in the text if they are used as math- 
ematical symbols. Punctuation marks are used at the end of equations as if they 
appeared directly in the text. 

Theorem 1. Theorems, lemmas, etc. are to be numbered consecutively in the 
paper. Use double spacing before and after theorems, lemmas, etc. 

Proof. Proofs should end with □ 
4. Illustrations and Photographs 

Figures are to be inserted in the text nearest their first reference. Where possible, 
figures should be embedded in electronic form from scanned images or a computer- 
based drawing package. Otherwise, original india ink drawings of glossy prints are 
preferred. Please send one set of originals with copies. If the author requires the 
publisher to reduce the figures, ensure that the figures (including letterings and 
numbers) are large enough to be clearly seen after reduction. If photographs are to 
be used, only black and white ones are acceptable. 



Fig. 1. A schematic illustration of dissociative recombination. The direct mechanism, 4mJ is 
initiated when the molecular ion Sl captures an electron with kinetic energy. 

Figures are to be sequentially numbered in Arabic numerals. The caption must 
be placed below the figure. Typeset in 8 pt Times roman with baselineskip of 10 pt. 
Use double spacing between a caption and the text that follows immediately. 

Previously published material must be accompanied by written permission from 
the author and publisher. 
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5. Tables 



Tables should be inserted in the text as close to the point of reference as possible. 
Some space should be left above and below the table. 
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Table 1. Comparison of acoustic for frequencies for piston-cylinder 
problem. 



Piston mass 


Analytical frequency 
(Rad/s) 


TRIA6-S1 model 
(Rad/s) 


% Error 


1.0 


281.0 


280.81 


0.07 


0.1 


876.0 


875.74 


0.03 


0.01 


2441.0 


2441.0 


0.0 


0.001 


4130.0 


4129.3 


0.16 



Note: Table notes 
Table footnote A 
''Table footnote B 



Tables should be numbered sequentially in the text in Arabic numerals. Captions 
are to be centralized above the tables. Typeset tables and captions in 8 pt Times 
roman with baselineskip of 10 pt. 

If tables need to extend over to a second page, the continuation of the table 
should be preceded by a caption, e.g. 11 Table 2. (Continued)" 
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References in the text are to be numbered consecutively in Arabic numerals, in the 
order of first appearance. They are to be typed in superscripts after punctuation 
marks, e.g. ". . . in the statement.^' 

References are to be listed in the order cited in the text. Use the style shown in 
the following examples. For journal names, use the standard abbreviations. Typeset 
references in 9 pt Times roman. 

7. Footnotes 

Footnotes should be numbered sequentially in superscript lowercase roman letters^ 
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This section should come before the References. Funding information may also be 
included here. 

Appendix A. Appendices 

Appendices should be used only when absolutely necessary. They should come after 
the References. If there is more than one appendix, number them alphabetically. 
Number displayed equations occurring in the Appendix in this way, e.g. (|A.1[> . 
(A.2), etc. 



a Footnotes should be typeset in 8 pt Times roman at the bottom of the page. 
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fj,(n,t) =^l(di <t,N(di) =n) l(N(cr) = n)da . (A.l) 

i=l Ja=0 
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